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Abstract: We present a new formula for all single trace tree amplitudes in four di¬ 
mensional super Yang-Mills coupled to Einstein supergravity. Like the Cachazo-He-Yuan 
formula, our expression is supported on solutions of the scattering equations, but with 
momenta written in terms of spinor helicity variables. Supersymmetry and parity are 
both manifest. In the pure gravity and pure Yang-Mills sectors, it reduces to the known 
twistor-string formulae. We show that the formula behaves correctly under factorization 
and sketch how these amplitudes may be obtained from a four-dimensional (ambi)twistor 
string. 



1 Introduction 


In any space-time dimension, the scattering equations underpin the classical S-matrix of a 
wide variety of massless field theories by constraining external kinematic data in terms of 
marked points on an auxiliary Riemann sphere, S [1]. Geometrically, these constraints can 
be compactly summarized as the requirement that a meromorphic quadratic differential 
P'^{z) vanishes globally on S [2], Recently, the scattering equations have received consid¬ 
erable attention for their central role in the Cachazo-He-Yuan (CHY) expressions for the 
tree-level S-matrices of various massless bosonic field theories [3, 4], although they were 
first discovered in the context of high-energy string scattering [5-8]. 

Even before the advent of the CHY formulae, the importance of the scattering equa¬ 
tions was realized for field theory in four space-time dimensions [9]. In four dimensions, 
the spinor helicity formalism allows us to solve the on-shell condition for momentum by 
writing it as the product of two Weyl spinors. In addition, the simplicity of on-shell super¬ 
space makes it possible to account for arbitrary helicity external states in super-Yang-Mills 
theory and supergravity. Given the existence of such expressions for the pure gravity and 
gauge theory sectors [10, 11], it is natural to ask if there is a generalization to the tree-level 
S-matrix of supergravity coupled to super-Yang-Mills theory (sEYM) in four dimensions 
(we consider the case where the spin one particles of the gravity multiplet remain un¬ 
gauged [12]). In this paper, we propose such a formula for the coupling of a single trace of 
gluons to gravity. It is compact, written in terms of integrals over the moduli space of a 
punctured Riemann sphere S, and supported on solutions to the scattering equations. Eur- 
thermore, it easily incorporates supersymmetry thanks to the simplicity of four-dimensional 
on-shell superspace, is manifestly parity invariant, produces the three-point amplitudes of 
the EYM action, and factorizes appropriately. 

The new formula suggests a remarkable structure in four-dimensions that would be 
completely obscured by a naive replacement ki ■ kj —)■ [ij]{ij) in the Pfaffians of the CHY 
formula. Remarkably, as in [11] the [, ] and (, ) factors completely decouple, appearing 
in separate determinants. The previously known expressions for amplitudes in pure gauge 
theory or supergravity can be derived from the sphere correlation functions of certain 
worldsheet theories [13-16] (collectively referred to as ‘twistor string theories’), and our 
formula also appears to have a worldsheet origin. 

2 The Formula 

Before describing our formula, we first point out several ab initio constraints which any 
purported formula for sEYM must obey. The first of these is rather obvious: in the pure 
gauge theory or supergravity sectors, it must reduce to known expressions for these tree- 
level S-matrices. Such expressions are provided by the Roiban-Spradlin-Volovich-Witten 
(RSVW) formula for super-Yang-Mills theory [9, 10], and the Cachazo-Skinner (CS) for¬ 
mula for supergravity [11]. 

A second constraint comes by considering the gravitational coupling constant k ~ 
V^Gn- a tree-level amplitude in EYM with n external gravitons and r colour traces must 
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be proportional to In particular, in the purely gravitational sector this is the 

usual K^~‘^ factor associated with (the Euler character of) a gravitational tree graph, while 
if there are no external gravitons and only a single colour trace, then we find as expected 
for the conformally invariant Yang-Mills answer. In a connected tree, different colour traces 
interact by exchanging gravitons, with a gravitational coupling at each end. 

In [11] it was explained that, when written in terms of a worldsheet model, these powers 
of K must be balanced by the same number of powers of [, ] or (, ) brackets (interpreted as 
infinity twistors, representing the breaking of conformal invariance). Parity transformations 
exchange [, ] and (, ), and so we learn that in four-dimensional EYM with n± external 
gravitons of helicity ±2 and r colour traces, we should obtain 


#(, ) = n_-hr - 1, #[, ] = n+-hr - 1. (2.1) 


on the worldsheet. 

A hnal constraint is given by a rather curious observation about the tree amplitudes 
of sEYM theory in four dimensions: every colour trace of external gluons must contain at 
least one gluon of each helicity.^ In particular, amplitudes with all negative helicity gluons 
and arbitrarily many positive helicity gravitons must vanish, despite the fact that this is 
far from obvious by standard MHV counting. To prove this, first note that each three-point 
interaction of the sEYM Lagrangian must contain at least one gluon of each helicity, or 
no gluons at all. This clear for the pure gauge theory interactions, and all three-point 
interactions with two gluons and one graviton arise from the term 

gP"" d[f,Ap] , 

in the Einstein-Yang-Mills action. Spinor helicity variables make it easy to see that any 
three-point function coming from this term with two gluons of the same helicity vanishes. 

Now, consider any tree diagram in sEYM, and select all of the gluons in the diagram 
belonging to a particular colour trace. This identifies a unique tree sub-graph associated 
with the chosen trace; we may assume that all its internal edges are gluon propagators, since 
a graviton propagator would lead to a double trace or an external graviton contribution 
which can be amputated. The total number of gluons (both internal and external) of each 
helicity appearing in this sub-graph is equal to the number of internal edges (helicity- 
labelled gluon propagators) plus the number of external gluons of that helicity. The Euler 
characteristic of the tree sub-graph, combined with the fact that the number of gluons of 
each helicity at every vertex is positive, then implies the claim: the number of external 
gluons of each helicity in the trace must be strictly positive. 

We now turn to the expression for the tree amplitudes, beginning with the non- 
supersymmetric case. In this paper, we will only consider the case with r < 1. External 
particles are specified by two spinors, Xia and Xia for particle i, as well as a helicity label. 
We divide the external particles of each amplitude into sets of gluons and gravitons of 

^Properties of EYM closely related to this were observed in [17]. 
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positive or negative helicity; positive (negative) helicity gravitons form the set h (h), and 
positive (negative) helicity gluons form the set g (g). 

Let us first present the formula, and then explain its various ingredients. The tree-level 
amplitude with a single colour trace is given by an integral 


fdet'$det'4 ^ 

“ yv„iGL(2,c) 


dti 


\2h\ 




dik 


fcsiiug 


\2h\-l 


5HXk- 


( 2 . 2 ) 


where h denotes the helicity of the given particle (i.e., h = ±1,2). The {zi^Zk} label 
marked points on an abstract Riemann sphere S = CP^ in an inhomogeneous coordinate, 
z, while the complex scaling parameters {ti,tk} carry conformal weight, taking values in 
expressions X{zi),X{zk) appearing in (2.2) are given by 


K^i) = h>^kS{i,k), Hzk) = ^ tjXjS{k,j), (2.3) 

fcehug j'ehug 

where 

Vdz dy 

S{z,y) := -, (2.4) 

z-y 

is the Szegd kernel at genus zero (free fermion propagator). Hence, the parameters tj, tk 
carry opposite weight with respect to little group scalings. 

The three main ingredients in (2.2) are the insertions det^^^, det^4>, and PT. $ is a 
(|h| +1) X (|h| ± 1) symmetric matrix with rows and columns corresponding to each positive 
helicity graviton and colour trace in the amplitude. Its entries are given by; 


^ij = titj[ij]S{i,j) fori/j, S{i,m), (2.5) 

meg 


^ij ‘hjg , <hgg — <hgj , 

jeh\{i} ieh 

where i,j G h, and g labels the row/column for the gluon trace. It is easy to see that <I> has 
co-rank one, with kernel spanned by the vector (1,1,... , 1), so its determinant vanishes. 
The operation det^ corresponds to removing any choice of row and column from a matrix 
and then taking its determinant. det^$ is easily seen to be independent of the choice of 
row and column removed, and is a canonically-dehned non-vanishing object. 

<I> is the parity conjugate of <I>. It is a (|h| ± 1) x (|h| ± 1) matrix, with entries 

^ki = tkii{kl)S{k,l) fork^l, ^kg = ikim {km) S{k,m), (2.6) 

meg 

^kk = - ^ ^kl- ^kg, ^gg = - ^ ^gk, 

ieh\{k} fceh 

where k,l G h. This matrix also has co-rank one with kernel spanned by (1,1,..., 1). 
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Finally, PT denotes a ‘generalized’ Parke-Taylor factor corresponding to the colour 
trace. If the gluon trace contains n total gluons (of all helicities), then this Parke-Taylor 
factor is 

n 

PT:= tr(TMi)...TM™.))]j5(^^(.),^^(.^,)), (2.7) 

O'^Snf’^n ^ 1 

where are the generators of the gauge group. 

Although takes the form of an integral expression, the delta functions in (2.2) 

saturate all of these integrals in addition to providing overall momentum conservation. So 
in reality, all integrals in (2.2) are performed algebraically against delta functions. These 
delta functions are a refinement of the usual scattering equations; they imply ki ■ P{zi) = 0 
for Paa{z) = Xaiz)Xa{z) [9]. Note that in the special cases where h = h = 0org = g = 0it 
is equivalent to the RSVW or CS formulae respectively, presented in the guise of [16]. The 
reduced determinants and definitions (2.5) and (2.6) ensure that is consistent with 

the counting of (2.1) and that the colour trace contains at least one gluon of each helicity. 

Supersymmetry can be incorporated straightforwardly, unlike in the CHY formulae, 
due to the simplicity of on-shell superspace in four dimensions. Extending EYM to AA < 3 
supersymmetry,^ on-shell scattering states are specified by the usual two Weyl spinors 
AQ,,Aa as well as Grassmann parameters for the supermomentum, or fiA, for A = 
1,... ,J\f. Remarkably, our formula accommodates this supersymmetry with the inclusion 
of a single exponential factor: 


M^’~ = 

h,h 


det'd> det'i 
vol GL(2, C) 


/ 


PT exp 




^ tiikfiAiTik S{k,i) 


iehUg 

VfcShUg 


/ 


X 




iehUg 


fcehUg 


\2h\-l 


( 2 . 8 ) 


As usual, amplitudes for individual helicity components of the supermultiplets are read 
off by expanding the exponential and extracting those terms of appropriate degree in the 
Grassmann variables. Note that both (2.2) and (2.8) are manifestly parity symmetric. 


3 Justification 

In this section, we show that (2.8) factorizes appropriately and produces the correct three- 
point amplitudes. 

3.1 Three-point amplitudes 

In the helicity-based framework, all three-point amplitudes of the EYM Lagrangian are 
classified as MHV or MHV, depending on whether they have one or two positive helicity 

^We believe that the formula is also correct for A/” = 4, provided one chooses an appropriate representa¬ 
tion for the external gluons. 
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legs, respectively. Since the formula (2.2) is parity-symmetric, we only check the MHV 
three-point amplitudes explicitly. For Yang-Mills gauge group SU(A^) there are only three 
such potential amplitudes: 

( 12)6 _ _ ^ 12)3 ^ _ (2 3 )^ 


hi h2 hs 


(2 3 ) 2 ( 31 ) 2 ’ ■ (2 3 )( 31 ) ’ - ^ 12)2 ’ 

with overall momentum conservation and any colour trace stripped off. Note that one can 
also write down an expression with the homogeneity of a gig 2 h 3 amplitude: 

- - (12)^ 
glg2h3 ~ (2 3)2(31)2 ’ 

but this does not occur in EYM (as follows from dimensional analysis). 

Since (2.2) reduces to known expressions in the all-gluon or all-graviton sectors, it 
is obvious that it reproduces the first two amplitudes in (3.1). So the only non-trivial 
calculation is to make sure that (2.2) produces the third amplitude in (3.1). In this helicity 
configuration, our formula becomes 

dt2 dh 


gi g2 h3 


(3.1) 


tr(T"iT"2 


) J —(2 3) 5^(Ai — ^2^2 — ts^s) S ‘^(\2 + i 2 ^lJ ) 


where the positions of zi,Z 2 , Z 3 have been fixed with the SL(2, C) freedom, the C* freedom 
has been used to set = 1, and t 2,3 have been rescaled. The final two integrations can be 
done against the delta functions in a straightforward manner, leaving 

/ 3 \ 


tr(rir2) 


^ 2 = 1 


(2 3)^ 
0^ 


as required. 


3.2 Factorization 

For the (non-supersymmetric) expression (2.2) to factorize correctly, we must show that in 
the limit where a subset L of the external momenta go on-shell 

lirn M({Ak^k,hk}) = 

(ZligL AiAi)2^0 

f d2A d2A 

/i=± 

where the sum is over possible helicities flowing through the factorization channel, the 
integral is over the on-shell phase space of the intermediate state, and i? is the compliment 
of L. In many respects, this calculation follows similar lines to those for the RSVW and 
CS formulae [18, 19]. 

In terms of the Riemann surface S underlying (2.2), the factorization limit should 
correspond to a degeneration of S into two Riemann spheres and joined at a node. 
Locally, we can model this degeneration in terms of inhomogeneous coordinates by 

(z - Za)(z - Zb) = q , (3.3) 
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where in the g —)• 0 limit the node is located at Za S Sl and One advantage of 

this local model is that the behaviour of the propagator (2.4) is particularly simple near 
the degenerate limit: 




S{i,j) if Zi,Zj G Si or Zi,Zj G Tjr 
-^^—S{i,a)S{b,j) +0{q) if 2:* G Zj G T.r 


(3.4) 


This follows from the universal behaviour of Szegd kernels on degenerate Riemann sur¬ 
faces [20]. 

Since the scaling parameters ti, carry conformal weight, their behaviour in the 
degeneration limit is non-trivial. The local model (3.3) dictates that a section of T-^ will 
scale as in the g —)• 0 limit, depending on which side of the degeneration the section 

is located.'^ The choice of which of or S/j is associated with the ‘-I-’ scaling must be 
summed over; below we see that this choice is associated with the helicity (positive or 
negative) of an intermediate particle flowing through the factorization channel. 

Homogeneity of the measure on the moduli space combined with little group scaling 
requires that the ti and tk parameters scale oppositely in the degeneration limit. The 
behaviour of the scaling parameters is thus given by: 

ti for Zi £T,l ~ j 4 for Zk G Sl 

g"Fi/4 f. for Zi G'Er ^ | ik for Zk G Y^r 



with the choice of upper or lower sign to be summed over. 

It is natural to work in a formalism where the only objects carrying conformal weight on 
S are these parameters and the Szego kernels. This is accomplished by defining parameters 
t*, 4 (valued in respectively) via a single insertion of 


1 = 


d4d4 ^ 
vol C* 



y/dZa dZb 


(3.6) 


in the amplitude (2.2). This allows us to re-write the behaviour of the propagator (2.4) in 
the attractive form: 




S{i,j) if Zi,Zj G Yl or Zi,Zj G Yr 
y/qUi* S{i,a) S{b,j) + 0{q) if Zi G Yr, zj G Yr 


(3.7) 


A priori, the t*, 4 are just convenient dummy variables, but they will eventually become 
associated with an intermediate on-shell particle flowing through the degeneration. 

It is also convenient to take a factor of or from each graviton wave function 
and incorporate it into det^<I> or det^<I>, respectively. Additionally, we can choose a single 

®To be precise, the local model should read {zl — Za){zR — Zb) = q, where zljZr are appropriately chosen 
coordinates on We keep this choice of local coordinates implicit in what follows to streamline 

notation. 

'^In actuality, (3.3) only restricts a section of scale as on El and qT(2-c«)/4 

consider the symmetric case q = 1 for simplicity only. 
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gluon of each helicity, say r £ g, s € g, and divide the trace row and column in each of 
$ by the associated scaling parameter, tr, ig- 

The result of these (trivial) manipulations is a transformation 


det'$ det'$ ^ 


iSh 


nf 

k&h ^ 


Pg det'$ det'l> 


ti 


n 

feeh 


dtfc 

tk 


(3.8) 


inside (2.2). Here, we abuse notation by writing and for the rescaled matrices with 
entries 


^ij = [ij]S{i,j) for i^j, [im] S{i,m) 


meg 




(3.9) 


‘^- = - E = - E r’ 

Ln 67 Lr 

jeh\{i} iGh ^ 


and likewise for <f>. These rescaled matrices have the (important) virtue of behaving nicely 
in the factorization limit. The definition of the reduced determinants is adapted appropri¬ 
ately for this rescaling: 


det'^> = 

ti tj 

det'<l> = 

tkk 


\^l\ 



where denotes the determinant of with row a and column b removed, etc. 

Note that the formula does not depend on which representative gluons we choose 
for rescaling the matrices. We always have at least one gluon of each helicity and the 
fundamental properties of the reduced determinant make it obvious that the final answer 
does not depend on the choice we make. Additionally, the rescaling of the matrix entries 
is reflected in a rescaling of the null vector; for instance, the kernel of $ is now spanned by 
(H,... ,f|h|)ir)) which scales according to (3.5) in the degeneration limit.^ 

Clearly, there are two different ways in which the formula (2.2) can factorize as g —)• 0: 
the degeneration may or may not split the colour trace. We will show that in the former 
case the intermediate state corresponds to a gluon, while in the latter it will be a graviton. 

Let us begin by considering the case which does not disturb the colour structure. The 
behaviour of (2.2) in the g —)■ 0 limit can be easily expressed by arranging the matrices 
$ and $ in a judicious manner. Arrange the entries of both matrices into blocks, so that 
the upper-left block corresponds to entries with both indices on Ej;, and the bottom-right 
block corresponds to entries with both indices on T,r: 


y^RL ^RR j ’ \^RL ^RR J 

^After the various rescalings, the reduced determinants both scale like O(y^) in the g —>■ 0 limit; to 
make this manifest in subsequent computations, we always choose to remove a row/column corresponding 
to an entry of the null vector which does not tend to zero under the degeneration. 
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Without loss of generality, we assume all gluons to remain on In the q ^ 0 limit, (3.7) 
ensures that the off-diagonal blocks of the rescaled matrices vanish at O(y^), so we must 
focus on what is happening in the diagonal blocks. Consider clearly the entries {^LL)ij 
and {^RR)ij are unchanged on and T,r, respectively, as g —)• 0. However, the diagonal 
entries behave as 

{^LL)ii = - ^ '^^[ij]S{i,a)S{b,j)+0{q), 

Lj Li Li 

i^RR)ii = - ^i^RRh S{i,b)S{a,j) + 0{q ^/‘^), 

jehi{\{i} * ieh^ug * 

where we choose the upper sign in (3.5) for concreteness. As the worldsheet factorizes, the 
particles h/j on T^r form an effective particle insertion on S/, which appears on the diagonal 
entries of while all entries of ^rr only encode the particles of hjj. The interpretation 
is clearly that the node becomes the insertion point of a new particle with positive/negative 
helicity on (The opposite configuration follows by taking the lower sign in (3.5), and 

the hnal result contains the sum over both choices.) 

To make this manifest, introduce a new spinor-helicity variable A* by inserting 

1 = J - UX{zb)^ (3.10) 

into the amplitude (2.2) close to the degeneration, with 

H^b) ■■= tj Aj S{b,j). (3.11) 

iShfl 

On the support of this delta function, the diagonal entries of d> can be rewritten as 
{^LL)ii = - ^ ■^{^LL)ij - ■f{^LL)ig + 0 {q) , 

i6h*\{i} * 

{^RR)ii = - ^ + 0 ( 9 ) • 

where = h^ U {*}. This is precisely the form required for the expected channel. 

With our choices corresponding to a positive helicity intermediate state on Tir, it is 
natural to remove a row and column corresponding to a graviton on 'Er when computing 
the original det^d*. As g —)• 0, one finds 

det'd* = y/q det'^L det'^R + 0{q), (3-12) 

where is the (|h)^|-|-l)x(|h2|-|-l) matrix appropriate for with the positive helicity 
intermediate state * at S S/,, and ^r is the |hij| x |h/j| matrix appropriate for S/j. 

The story for <I> proceeds in a similar fashion. By introducing 

1 = y d^A* (5^(A* - t*A(za)) , (3.13) 
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in (2.2), with 


(3.14) 


-^(^a) := ^ ikXkS{a,k), 

fcehiUg 

the diagonal entries of ^ll and ^rr become 

{^LL)kk = - y2 ^ ^kg + 0{q ), 

lehL\{k} 

{^RR)kk = - ^ -Y^ki + 0{q), 

i&h*^\{k} ^ 

where h|j = h/j U {*}. With this corresponding to a negative helicity intermediate state on 
Sjj, we take det^'h by eliminating a row and column for a graviton on T,r, to obtain the 
desired factorization: 

det'4> = det'^L det'^j? + 0{q). (3.15) 

This shows that the reduced determinants factorize correctly and yield a factor of qt'^i'^, 
while also introducing the wave functions for an intermediate on-shell state. 

Next we examine the behaviour of the external wave functions. On the support of 
(3.10), (3.13), the arguments of the delta functions for the external states involve 


tiX{zi) = 1 

f T.k&hLUetitkXkS{i,k)+ 0{q) 

[ Xyfceh|j ^*4 Afc S{i, k) 

for Zi G Si 
for Zi e'SR 

(3.16) 

4 A(zfc) = 

i X/ieh^ug ^*4 Aj >S'(A:, i) 

for Zk G Sl 

(3.17) 

\ SiehflU ^*4 Xi S{k, i) -|- 0{q) 

for Zk G T,r 

in the degenerate limit g —)• 0. 




As g —>• 0 a scaling argument based on the local model (3.3) dictates how the measure 
on the moduli space of S factorizes into measures on the moduli spaces of and Tir (cf. 
[21, 22]). Combining all of our observations up to this point, the formula (2.2) looks like 


S 1 - 


d(7d^A*d^A* dzadzb 

7 (volC*)2 (vo1SL(2,C))2 

X det'<hL det'$L det'<l>i^ det'^i? PT 


-t^X{Za)) 5"^ {Xi - tiX{Zi)) 


dife r 2 


^ 6 ^{X,-U 


iehLUgi 

dt 


- t 


\2h\-l 


fcehiUgi k 

dii 


n n r\2h\-l 


j&iiR 


lehR 


(3.18) 


up to 0 {q), where we have reverted to the original, unrescaled matrices $,$. 

The delta functions for A* and A* are naturally incorporated into products over Oieh* ug 
and Ozeii* gravitons. Finally, trading the delta function of (3.6) for an additional 

R 
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(vol C*) ^, we are left with 


where 


f dq d^A* d^A* 

J q vol C* 


{MtMJ, + 0 (q)) , 


Mt = 


det'<f>Ldet'^L 
vol GL(2,C) 


pT n 

ieh* Ug^i 


UHzi)) J2h|-l' 

fcehiUg k 


(3.19) 



and similarly for A4^. 

All that remains to show is that extracting the residue at g = 0 corresponds to setting 
the momentum flowing from to T,r on-shell. To this end, notice that the various delta 
functions in (3.19) imply that 


AjAj -|- AfcAfc — A*A*-|-0(g). 

iehiUg fcehiUg 

In addition to enforcing momentum conservation on each side of the cut, this reveals the 
intermediate particle as manifestly on-shell. Thus, taking the residue at g = 0 results in 
the h = + graviton exchange term in the factorization expression (3.2). The h = — term 
is obtained in a similar fashion, by choosing the other sign in (3.5). 

The second possible degeneration, which disturbs the colour structure of the ampli¬ 
tude via a gluon exchange, follows in much the same way, so we will be more brief in its 
description. In this case the external gluons are split between and S/j in the g —>• 0 
limit. The behaviour of the matrices ^ and under the degeneration is the same as before; 
it is now convenient to eliminate the row and column corresponding to the trace in both 
det'<I> and det'$. 

Once again, after inserting the delta functions (3.10), (3.13), the reduced determinants 
factorize as (3.12), (3.15) with the intermediate particle being incorporated into = 
gL U {*} on Si and = g/jU {*} on S/j. However, the factor of qt^t^ is entirely absorbed 
by the scaling of the gluon representatives appearing in (3.8). 

From (3.7) it is easy to see that the Parke-Taylor factor behaves as 

PT = q tlP, PTi. PT^, + 0{q ^/^), (3.20) 


where PT-y£ denotes the Parke-Taylor factor for those gluons located on S^ with the point 
Za inserted in the cyclic ordering precisely where the original colour trace is broken by the 
degeneration. The measures for appearing in (3.6) are automatically appropriate for 
gluons, so we obtain the correct products over delta functions for factorization. Taking the 
residue of the resulting dq/q measure corresponds to the h = + gluon exchange term in 
the factorization expression (3.2). 

Factorization for the supersymmetric amplitude (2.8) is established by looking at how 
the exponential encoding supersymmetry behaves in the g —>• 0 limit. For concreteness we 
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pick the case of (3.5) where ti ~ on and the other case follows analogously, as 
usual. To begin notice that the exponent splits as 


Uikfji-'nk S{k,i) + ^ tiikfii-'nkS{k,i) + ^ tiiktd* fji-r]k S{k,a)S{b,i) + 0{q). 

iehiUgi iehijUgH iehflUgH 

fceiliUgL fcehflUgH feeh^ugi 

(3.21) 

The first two terms are clearly appropriate for and respectively, while the third 
term accounts for the new intermediate particle. In fact, a simple calculation reveals that 
the exponential of this third term can be written 


/ 




exp 


tiikUi*fii-r]kS{k,a)S{b,i) 


. 2eh_RUg_R 
\/i;Gh_LUgL 


= y d-^? 7 *d-^? 7 * e’'*'^* exp | ^ tJkV*-VkS{k,a)+ ^ S{b,i) 


\kehLUgL 


i&hnUgu 


on the support of the same delta functions used in the purely bosonic calculation. 

In summary, the exponential encoding supersymmetry factorizes, with the arguments 
on Si, Sij becoming 


tdk fii-T]kS{k,i) and ^ tdk fji ' Vk S{k,i), 

ieh^ugi, iehflUgfl 

A:6hz,Ugz, fc6h|jUgfl 


respectively (without altering the simple pole in q). Combined with our previous argu¬ 
ments, this leads to the appropriate on-shell superspace measure 


dg d^l^A, d^l^A, 
q vol C* 


3'?*-^* (MjM 


R 


+ Oiq)) 


(3.22) 


So factorization of the super-amplitude follows straightforwardly from factorization of the 
bosonic formula. 


4 Conclusions 

We have presented a new formula for all single trace tree amplitudes in supersymmetric 
Einstein Yang-Mills theory in four dimensions, written in terms of on-shell superspace. 
The formula was shown to reproduce the correct three-point amplitudes and to factorize 
appropriately in both gravitational and coloured channels. The considerations at the be¬ 
ginning of Section 2 also place strong constraints on the form of general multitrace EYM 
amplitudes. It would be interesting to investigate these further. 

In the purely gravitational sector, the formula reduces to the representation of am¬ 
plitudes given in [11, 16]. These formulae are known to be the output of a twistor-string 
theory for maximal supergravity in d = 4 [15] and it is natural to wonder whether there is 
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a modification of this theory that describes sEYM. In this regard, we note that the cou¬ 
pling of the gluons to the gravitons in d>, together with the Parke-Taylor factors, may 
be generated by inserting operators 

, and iv (p-^ Oa D-^ O. 


Here D = d + A{Z) + A{W) and 


r... dA] 


w — 

’ dZ 

+ 


- K 

dZdZ^ 




dA 

'm 


/ _ \ 

+ pP^dWdWK / ’ 


where A and A are gluon wavefunctions. It seems likely that the RSVW formula for tree 
amplitudes in sYM are best interpreted as the single trace sector of a twistor string for 
sEYM. Such a theory would also enable the use of worldsheet factorization arguments to 
streamline the calculations above [23]. 
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